PHYSICAL REVIEW D 68, 123521 (2003) 



Observational constraints on Chaplygin quartessence: Background results 



o 
o 

(N 



Martin Maklcr, Sergio Quinet de Oliveira, and Ioav Waga 

Universidade Federal do Rio de Janeiro, Instituto de Fisica, 
Caixa Postal 68528, CEP 21941-972 Rio de Janeiro, RJ, Brazil 
(Dated: Received 27 June 2003; published December 2003) 

We derive the constraints set by several experiments on the quartessence Chaplygin model (QCM). 
In this scenario, a single fluid component drives the Universe from a nonrelativistic matter-dominated 
phase to an accelerated expansion phase behaving, first, like dark matter and in a more recent epoch 
like dark energy. We consider current data from SNIa experiments, statistics of gravitational lensing, 
FR lib radio galaxies, and x-ray gas mass fraction in galaxy clusters. We investigate the constraints 
from this data set on flat Chaplygin quartessence cosmologies. The observables considered here 
are dependent essentially on the background geometry, and not on the specific form of the QCM 
fluctuations. We obtain the confidence region on the two parameters of the model from a combined 
analysis of all the above tests. We find that the best-fit occurs close to the ACDM limit (a = 0). 
The standard Chaplygin quartessence (a = 1) is also allowed by the data, but only at the ~ 2a 
level. 
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I. INTRODUCTION 

Over the past decade, a cosmological model consis- 
tent with most astrophysical data available to date has 
emerged: a flat Universe whose evolution is dominated 
by a repulsive cosmological term [or dark energy (DE)] 
and pressureless cold dark matter (CDM). In addition to 
observational evidence for this model, there are also the- 
oretical motivations for it: for example, inflation theory, 
which generates a nearly flat space geometry and scale 
invariant primordial perturbations. Another example re- 
gards the nature of the two dark components. Particles 
predicted by extensions of the standard model of interac- 
tions, such as the lightest supersymetric (SUSY) particles 
or the axion Q are, for instance, natural candidates 
for CDM. On the DE side, a slowly rolling scalar field 
has been known to produce accelerated expansion since 
the proposal of inflation. It is thus a well motivated can- 
didate for the cosmological term ||. In fact, the sim- 
plest and most popular candidate as the driving force for 
the accelerated expansion is the cosmological constant 
A. However, its tiny value inferred from present obser- 
vations creates a major puzzle. If A is to be considered 
as the sum of the vacuum energies of all fields, it is hard 
to understand how they would cancel to one part in 10 55 
or 10 122 (for a review on the cosmological constant and 
DE, see ||). 

In the standard cosmological model, dark matter and 
dark energy are necessary to account for two seemingly 
independent phenomena — clustering of matter and ac- 
celerated expansion. As we pointed out, there are some 
theoretical hints about the nature of these two compo- 
nents, but, in fact, at present, there is no conclusive ob- 
servational evidence that these phenomena are produced 
by distinct components (see a first attempt to address 
this question in Ref. Therefore, instead of using 

a theoretical bias towards dark energy plus CDM, one 
may choose an alternative point of view and search for 



phenomenological models that could still be consistent 
with current observations, motivating theoretical investi- 
gations a posteriori. In addition to being potential candi- 
dates for the dark sector, these models would allow us to 
test the robustness of observational predictions regarding 
the determination of the cosmological model. 

From the point of view of simplicity, it would be inter- 
esting if instead of two unknown components we could 
have a single one, accounting for the phenomenology not 
associated with ordinary matter. In this model, a single 
component would be responsible for both clustering and 
accelerated expansion. Such a model is usually referred 
to as UDM (unifying dark matter energy) or — since there 
is only one dark component besides baryons, photons, 
and neutrinos — as quartessence 0. 

The possibility of a unified description of DE-DM has 
given rise to rather widespread interest recently The two 
major candidates for UDM are the quartessence Chap- 
lygin fluid (QCM) 0, |3 and a quartessence tachyonic 
field 0. Some confrontations of QCM against obser- 
vational data have been performed. It has been shown 
that, for a wide range of parameters, this model is con- 
sistent with a number of tests of the background metric 

Nevertheless, preliminary analyses of large-scale struc- 
ture and CMB data favor the ACDM limit of the QCM 
model |5lll5j. However, these perturbation analyses need 
further assumptions beyond the background Chaplygin 
equation of state. For instance, they assume that there 
are no viscous stresses in the perturbed fluid. In [T^ 
it was shown that, if entropy perturbations are allowed, 
instabilities and oscillations, present in the mass power 
spectrum in the adiabatic case, may be eliminated and, 
as a consequence, the parameter space is enlarged. 

The main goal of this paper is to set constraints on 
the Chaplygin quartessence model and check whether it 
is consistent with present cosmological data. We focus on 
observables that are dependent essentially on the back- 
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ground geometry. We perform a combined analysis of 
data, including SNIa experiments, gravitational lensing 
statistics, FRIIb radio galaxies, and x-ray gas mass frac- 
tion in galaxy clusters. Some of these tests were studied 
previously within the QCM setting. We review these re- 
sults with a careful treatment and discuss the outcome 
of a combined analysis of the data. 

The paper is organized as follows. In the next sec- 
tion, we discuss the phenomcnological motivation for 
quartessence and present its realization through a fluid 
with an exotic equation of state, focusing on the QCM. 
Constraints from SNIa experiments are discussed in Sec. 
Tu\ In Sec. IIVI we discuss the bounds on QCM param- 
eters imposed by the statistics of gravitational lenses. 
Limits from Fanaroff-Riley type lib radio galaxies are 
obtained in Sec. and from x-ray gas mass fraction in 
galaxy clusters in Sec. IVII We present a combined anal- 
ysis of the experiments and our concluding remarks in 

sec. ivra 



THE QUARTESSENCE CHAPLYGIN 
MODEL 



II. 



As discussed in the previous section, the Universe is 
believed to be dominated by two unknown components, 
generically denoted by dark energy and dark matter. At 
the cosmological level, the direct detection of each of 
these two components involves observations at different 
scales. Since it is not supposed to cluster at small scales, 
the effect of dark energy can only be detected over large 
distances, where the accelerated expansion is observed. 
On the other hand, the CDM is detected by its local clus- 
tering through the motion of visible matter or through 
the bending of light in gravitational lensing. 

Within the standard lore of general relativity and the 
Friedmann-Lemaitre-Robertson-Walker model, for the 
Universe to undergo acceleration its average density and 
pressure must obey 



( P + 3p) <0, 



(1) 



so that there is an effective gravitational repulsive effect. 
On the other hand, for the large-scale structures we see 
to have formed today, the dark matter has to be nonrel- 
ativistic, i.e., 



\p\ <p. 



(2) 



These two conditions are not in contradiction, since ob- 
servations in various scales probe different average densi- 
ties. For example, local motions in clusters of galaxies oc- 
cur in regions hundreds of times denser than the average 
density of the Universe. So, if the pressure is negative [as 
needed from Eq. and the ratio \p\ / p is a decreasing 
function of the density, the two conditions can be made 
compatible. In this picture, the quartessence would act 
as dark energy in very low-density regions and as dark 
matter in higher-density regions. 



A very simple equation of state that has the properties 
discussed above is the so-called generalized Chaplygin gas 

0,E2, 1,0, 



m 4(q+1 



(3) 



where M has the dimension of mass. Consider now the 
background geometry of the Universe, which can be de- 
termined by assuming a homogeneous fluid whose den- 
sity and pressure are given by the averaged values of the 
present clumpy matter distribution. In this case, the en- 
ergy conservation equation 



PCh 



(pch +Pch)3- 
a 



(4) 



can be easily solved. The energy density of the Chaplygin 
fluid will be given by 



PCh = PChO 



3(a+l) 



.4 



l/(a+l) 



(5) 



where a is the scale factor, ao is its present value, A = 
(M 4 /p c / l o)' Q!+1 \ and the overdot in Eq. |g} denotes the 
derivative with respect to cosmic time. The equation-of- 
state parameter and the adiabatic sound velocity, for the 
Chaplygin component, are given by 



wch(a) 



PCh 
PCh 



Aa 3( - a 



(1 - A) + Aa 3 (»+ 1 ) 



and 



PCh 



-sCh 



= - — = -awch(a) 

PCh 



(6) 



(7) 



In principle, the constant M can take any positive 
value, but we should have M ~ 10 -3 eV in order to 
have negative pressure at recent times 0. Note that, 
in QCM, we need A > 1/(3 — 3^60 ) to have cosmic 
acceleration starting before the present time (qo < 0). 
For instance, for Q^q = 0.04, the above condition im- 
poses the constraint A > 0.347. More generically, it is 
necessary that A > in order to have a positive cos- 
mological constant at late times. The adiabatic sound 
speed maximum value (which occurs in the regions where 
PCh — > — Pch) is given by ^fot. Therefore, the parameter 
a is restricted to the interval a < 1. The Chaplygin 
gas a = 1 is the extreme case, where the sound velocity 
can be nearly the speed of light. Another special case is 
a = 0, which gives the equation of state pch = —pchoA 
and pch = Pcho^4 + PCho(l — A)(a /a) 3 , and is, therefore, 
equivalent to a superposition of CDM and a cosmologi- 
cal constant. Note that, for a < —1, we have at early 
times wch ~ — 1 and at late times wch 0. Thus, since 
we are interested in the quartessence scenario, where at 
early times i«ch ~ and at late times Wch ~ — 1, we 
also impose a > — 1. For p^h to be well defined and/or 
positive at early times, we also impose A < 1. Therefore 
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in QCM, the parameter a is restricted to the interval 
— 1 < a < 1, and A is restricted to < A < 1. We can 
see that, in fact, in QCM Eq. ^ interpolates between 
dark matter and dark energy as the average energy den- 
sity of the Universe changes. That is, when a/ao <C 1, 
we have pch oc a -3 and the fluid behaves as CDM. For 
late times, a/ao S> 1, and we get pch = —pch = — M 4 = 
const as in the cosmological constant case. 

Equation (0) provides the simplest example of 
quartessence. It has naturally a single-dimensional con- 
stant, since it is given by a power law, and allows us to 
solve the background energy conservation equation ana- 
lytically. Besides, it has only two free parameters. How- 
ever, many other Ansatze for equations of state (EOS) 
satisfying these criteria can be found. 

So far we have introduced the QCM from a purely phe- 
nomenological point of view. The EOS in Eq. was 
chosen as a simple toy model that could allow dark en- 
ergy/dark matter unification. A completely independent 
point of view, motivated by brane dynamics, led to the 
same EOS. Kamenshchik et al. Q proposed EOS in Eq. 
©, with a — 1 (the standard Chaplygin case) as a dark 
energy candidate. Only after was it realized that such 
fluid could naturally lead to dark energy /dark matter 
unification |||. Some possible motivations for this sce- 
nario from the field theory point of view are discussed in 
Refs. 0, H, • The Chaplygin gas appears as an effec- 
tive fluid associated with the parametrization invariant 
Nambu-Goto c?-brane action in a (d + 1, 1) spacctimc, 
and it can also be derived from a Born- Infold Lagrangian 
dlHj- The generalized Chaply gin EOS can be obtained 
from a complex scalar field Lagrangian with appropriate 
potential . The relation between a Chaplygin- like 
gas and the tachyonic scalar field was investigated in [2(| • 

In the rest of this paper, we will consider only the case 
of the Chaplygin fluid as quartessence. We will place 
constraints on the QCM parameters from several cos- 
mological tests which are dependent essentially on the 
background geometry. We emphasize that the resulting 
constraints depend on our choice of priors. For exam- 
ple, in Chaplygin quintessence it is commonly assumed 
that S^cdm ~ 0.3, while here we assume Ocdm = 0. 
Therefore, the limits on the parameters a and A will be 
different. 

A fundamental quantity related to the observables con- 
sidered here is the distance-redshift relation, given by 



where Hq is the Hubble constant, z = ao/a — 1 is the 
redshift, and f^ch, fi& are the density parameters of the 
Chaplygin fluid and baryons, respectively. 

Notice that following the idea of unification, we will 
not include an additional dark matter component. Thus, 
in Eq. ((5J) only the baryonic matter scales as (1 + z) 3 . 
Just for a <C ao does the Chaplygin component scale as 
CDM. In this case, we have an effective matter density 
parameter 1 



O 



eff 



n Ch (i~A) 1/il+a) + n b 



(10) 



Notice that, after radiation domination, in this sce- 
nario the Universe is always quartessence-dominated, 
thus avoiding the "why now" problem. However, there 
is still some fine-tuning to determine whether the aver- 
age equation of state turns from "matter like" to "energy 
like" (i.c, when the quartessence pressure begins to play 
an important role |||). 

Because after the decoupling epoch the QCM has neg- 
ligible pressure, it will cluster as ordinary CDM. In this 
model, only the unclustcrcd part would lately have ap- 
preciable negative pressure. 

Since observations of anisotropics in the CMB indicate 
that the Universe is nearly flat, and since the inflation 
paradigm predicts a flat geometry, we restrict the fol- 
lowing discussions to the zero curvature case, such that 

^Ch + Ofc = 1. 

The baryon density and Hubble parameters can be de- 
termined independently of the quartessence model. The 
observed abundances of light elements together with pri- 
mordial nucleosynthesis give Clbh 2 = 0.0214 ± 0.0018 
||. The Hubble Space Telescope (HST) key 
project result is h = 0.72 ± 0.08 @, where H := 
lOO/i Km/sec/Mpc. These bounds on f2&/i 2 and h are 
in agreement with a fit of purely Wilkinson Microwave 
Anisotropy Probe (WMAP) and other Cosmic Microwave 
Background (CMB) data [24[- Thus, throughout most of 
this paper we fix the Hubble and the baryon density pa- 
rameters at h = 0.72 and fib = 0.041. With Qb and h 
being determined by independent measurements, in the 
following sections we will investigate the constraints on 
the parameters A and a from several experiments, as- 
suming the above-stated priors. 



■(z) 



[H (z')]' 1 dz' , 



(8) 



where H is the Hubble parameter. In the QCM case 
(neglecting the radiation component), we have 



H{z)=H ({f! fc (l + 2) E 



n 



Ch 



A+(l-A){l + z) 3(1+a) 



1/2 



(9) 



III. TYPE IA SUPERNOVAE EXPERIMENTS 

Supernovae constraints on the Chap lygin models have 
been analyzed by many authors |l0j.|6l ITll ITEj . using 
different priors and simplifications. Here we repeat the 



[1] This is why we have used the parametrization f2^ f := (1 — A) in 
Ref. that corresponds to the effective matter density of the 
generalized Chaplygin gas for flch — 1 an< f ct = 0. 
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analysis presented in |6j, but we shall not impose any 
prior on the parameter A. 

The luminosity distance of a light source is defined in 
such a way as to generalize to an expanding and curved 
space the inverse-square law of brightness valid in a static 
Euclidean space, 



/ L \V2 
dL=[^^) = {l + z)r(z). 



\AitF) 



(11) 



In Eq. L is the absolute luminosity, T is the mea- 

sured flux, and r(z) is given by Eq. (JSJ). 

For a source of absolute magnitude M, the apparent 
bolomctric magnitude m(z) can be expressed as 



m(z) = M + 51og\D L 



(12) 



where Dl = Dl(z, a, A) is the luminosity distance in 



units of H , and 



M =M-51og#o + 25 



(13) 



is the "zero point" magnitude (or Hubble intercept mag- 
nitude) . 

We follow the Bayesian approach of Drell, Loredo and 
Wasserman [2j| and consider the data of fit C of Pcrl- 
mutter et al. [2(| with 16 low-rcdshift and 38 high- 
rcdshift supernovae. In our analysis, we use the following 
marginal likelihood: 



C(a,A) 



sV2ir 
A?] 



-9/2 



(14) 



Here 



^ (-Slog^ + ^ + rngf) 2 

(«! A )= 2 



i=l 
38 

E 



-51ogL>L + v + to| 



a 



2 

high,?' 



(15) 



where 



v(a,A) 



' 16 

E 



51ogD L (z 4 ,a,^) - m§5 r 



low, 2 



+ E 



5logD L (zi,a, A) - mg 



eflf 



! = 1 



high,z 



\ i= l CT low,i l= i CT high,z/ 



low, 2 



5 loge 



(16) 



(17) 



(18) 






a 



FIG. 1: Constant confidence contours (68.3% and 95.4%) in 
the (a, A) plane of the QCM parameters allowed by SNIa 
data, as described in the text. 



and 



2 

"high,? 



+ 



5 loge 



(19) 



The quantities m^" 



-.eff 
l B i 



cr m efi and a z are given 



in Tables 1 and 2 of Perlmutter et al. |26j . 

The results of our analysis for the QCM, marginalizing 
over the intercept, are displayed in Fig. 1. In this figure, 
we show 68.3 and 95.4 confidence level contours in the 
(a, ^4) plane. As we had observed in 0, current SNIa 
data constrain A to the range 0.4 < A < 0.9, but do 
not strongly constrain the parameter a in the considered 
range. 



IV. GRAVITATIONAL LENSING STATISTICS 

The statistics of gravitational lensing [2]| is one of the 
most traditional and widely used methods to constrain 
cosmological models, especially those with a cosmologi- 
cal constant. Since the publication of the first works on 
lensing statistics [28|, |25j |3(J almost 15 years ago, several 
studies constraining ACDM models |3lL l32l l33l | and some 
of its variants [34l l35l l3r| appeared in the literature. Al- 
though not without controversy, before 1998 the general 
belief was that models with Q m o ~ 1 and f^Ao ~ are 
preferred by lensing calculations. These results started 
to be challenged, a few years ago, when the high red- 
shift supernovae observations, in combination with CMB 
measurements, consistently indicated a nearly flat, low- 
density, and accelerated Universe. In fact, lensing esti- 
mates were not in strong conflict with the supernovae 
observations. However, it cannot be said also that they 
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were in comfortable agreement. There was room in the 
parameter space for concordance, mainly if A is dynami- 
cal, but the best-fit regions for each one of the tests were 
not in good agreement. 

With few exceptions, in the past most lensing analyses 
were based on optically selected quasars. One important 
concern in lensing investigations with optical sources is 
extinction in the lensing galaxies. The presence of dust 
in the lens galaxy could make the quasar images fainter 
and, as a consequence, difficult to detect. Radio selected 
surveys are immune to extinction. Besides reducing the 
total number of multiply imaged systems in an optically 
selected sample, extinction also affects its angular sep- 
aration distribution and this is very relevant in lensing 
analyses. In most studies with optically selected sources, 
extinction has been neglected. The motivation for this 
assumption is justified by the fact that early-type galax- 
ies that dominate the lensing statistics are believed to 
have little dust at the present epoch. However, the pos- 
sibility that the existence of dust in early-type galaxies at 
high redshifts (z > 0.5) could reconcile lensing calcula- 
tions with models with high A was suggested 10 years ago 
by Fukugita and Peebles [31} . Malhotra, Rhoads, and 
Turner [3Sj explored this possibility further and found 
evidence for it. They estimated a mean extinction of 
Am = 2 ± 1 magnitudes. By observing that statistical 
lensing analyses based on optical and radio observations 
were not in accordance, Falco et al. j3i| suggested that 
they can be reconciled if the existence of dust in E/SO 
galaxies is considered. However, they estimated that 
the required mean extinction is considerably lower than 
that estimated in j^. They obtained Am = 0.6 ± 0.4 
mag. In a subsequent work, Falco et al. |40j concluded 
that a substantial correction for extinction is necessary 
in any cosmological estimate using the statistics of lensed 
quasars. Their directly measured extinction distribution 
was consistent with statistical estimates from compari- 
son of radio-selected and optically selected lens surveys 

a. 

Another important source of uncertainty in gravita- 
tional lensing investigations is the velocity dispersion of 
early- type galaxies (cr* e ^)- The probability of lensing, or 
the optical depth (r), depends on the fourth power of 

al e ^ and, hence, is very sensitive to this quantity. In 
|4lj . Kochanek advocated a high value for the velocity 
dispersion, namely ai^ ~ 225 Km/sec. Lensing analyses 

that use relatively high values of ai^ predict low val- 
ues for A. The reason for this is simple: high A values 
imply larger distances and, as a consequence, also larger 
probability for lensing. In order to have r in high-A mod- 
els in accordance with observations, one possibility is to 
"compensate" its increase due to A by a reduction in 
cr* e ^ . However, this is not so easy. The angular image 
separation, another observable in lensing statistics, de- 
pends on the second power of ai . Therefore, a too low 
value of cr; may not fit the image separation distribu- 
tion observed in optically selected lens surveys. In fact, 



the main support for a high cr* value comes from the 
angular image separation distribution observed in these 
surveys. The observed value in the sample used in |4l| is 
(A8) ~ 1.6 arcsec. The sample of radio selected multiply 
imaged sources used by Falco et al. |39l | in their compar- 
ison with the optically selected ones also has a similar 
high (A6>). 

The situation described above changed considerably 
with the completion of the Cosmic Lens All Sky Survey 
(CLASS) sample of radio sources Recently, Chae 

and Chae et al. 01 performed a statistical lensing analy- 
sis with CLASS data. Under a well-defined selection cri- 
terion, they selected a total of 8958 radio sources from the 
CLASS. From these, a total of 13 systems have multiple 
images, giving a lensing rate of ~ 1/689. An important 
difference between the new data and previous ones is that 
in this sample the mean angular image separation is con- 
siderably smaller, (A0) ~ 1.2 arcsec. By assuming a sin- 
gular isothermal ellipsoid as the lens model and a "steep" 
faint-end slope of early-type galaxies luminosity function 
(pS e ^ = —1), they determined the velocity dispersion to 

be ai^ ~ 198 Km/sec. This value is significantly smaller 
than the best-fit value obtained in from optically se- 
lected surveys. By assuming a flat ACDM model, their 
likelihood analysis gives fl m o = O-^-o'W at the 68% con- 
fidence level. They also obtained results without fixing 
the curvature but fixing the equation of state of the dark 
energy (w) to be w = — 1, and/or fixing the curvature 
to be zero and taking w constant but not necessarily 
equal to —1. Their results indicate that, although lensing 
statistics is less restrictive than the magnitude-rcdshift 
test with type la supernovae, the two tests are in very 
good agreement. 

Although lensing analyses with radio sources are more 
reliable, in this work we still use optically selected 
quasars. We leave for future work the inclusion of ra- 
dio sources in our analysis. 

In the following, we briefly outline our main assump- 
tions for the lensing analysis using highly luminous 
quasars (3^, ■ We consider data from the HST Snap- 
shot survey [498 highly luminous quasars (HLQ)], the 
Crampton survey (43 HLQ), the Ycc survey (37 HLQ), 
the ESO/Licge survey (61 HLQ), the HST GO obser- 
vations (17 HLQ), the CFA survey (102 HLQ), and 
the NOT survey (104 HLQ) gj|. We consider a to- 
tal of 862 (z > 1) highly luminous optical quasars plus 
five lenses. The lens galaxies are modeled as singular 
isothermal spheres (SIS), and we consider lensing only 
by early-type galaxies, since they dominate the lens pop- 
ulation. We assume a conserved comoving number den- 
sity of early- type galaxies |47j . n e = no(l + z) 3 , and 
a Schechter form |4o j for the early-type galaxy popu- 
lation, n = J °° n* (L/L*) a * exp (-L/L*) dL/L*, with 
= 0.64 ±0.19 /i 3 10- 2 Mpc" 3 and a* = -1.0 ±0.09 
[44I |49| . We assume that the luminosity satisfies the 
Faber- Jackson relation [Hof. L/L* = (cr/cr*) 7 , with 7 = 4 
and adopt er» = 198 Km/sec. We emphasize that, even 
if we consider velocity dispersion of early-type galaxies 
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FIG. 2: Confidence contours in the (a, A) plane from the 
statistics of optically selected gravitational lenses. 



as low as cr» ~ 190 Km/sec, extinction has to be con- 
sidered in order to reconcile the statistical lensing anal- 
ysis using optically selected sources with those based 
on CLASS radio-selected sources. Further, if we use 
n* = 0.99±0.05/i 3 lCT 2 Mpc- 3 and a„ = -0.54±0.02, as 
obtained by Madgwick et al. [5lJ, it would be necessary 
to consider a higher amount of extinction. Previous anal- 
yses with optically selected quasars that use high values 
for cr» and do not consider extinction are not consistent. 

For SIS, the total optical depth (t) can be expressed 
analytically, t{z s ) = F/30 [d A (l + z s )f (cH^ 1 )- 3 , 
where z s is the source redshift, d A = r(z)/(l + 
z) is its angular diameter distance, and F = 
167r 3 n*(c J ff,7 1 ) 3 (cr»/c) 4 r(l + a» +4/ 7 ) ~ 0.016 measures 
the effectiveness of the lens in producing multiple im- 
ages [23. We correct the optical depth for the effects of 
magnification bias and include the selection function due 
to finite angular resolution and dynamic range |33l l4lj| . 
We assume a mean extinction of Am = 0.6 mag; this 
makes the lensing statistics for optically selected quasars 
consistent with the results of [44J . 

In Fig. 2, we show contours of constant likelihood 
(95.4% and 68.3%) arising from lensing statistics for the 
quartessence Chaplygin model. It is clear from the figure 
that lensing statistics weakly constrains these models. 
Only models with A > 1 are excluded by the lensing 
data. 



V. FANAROFF- RILEY TYPE IIB RADIO 
GALAXIES 

In addition to gravitational lensing statistics and the 
SNIa magnitude- redshift test, another useful method to 
constrain cosmological parameters is the classical angular 



size-redshift test. In this paper, we shall be concerned 
with the Fanaroff-Riley type lib (FR lib) radio gal axy 

f version of this test as proposed in [53 (see also |5J, 
IBH l57j ) . This test consists of a comparison of two 
independent measures of the average size of the lobe-lobe 
separation of FR lib sources, namely the mean size (D) 
of the full population of radio galaxies at similar redshift 
and the source average (over its entire life) size , which 
is determined via a physical model that describes the 
evolution of the source. The basic idea is that (D) must 
track the value of such that the ratio i?* = (D) /D* 
is independent of redshift. It can be shown that i?„ oc 
r (2/3/3+3/7)^ wngr e r is the comoving distance, and (3 is 
a parameter to be determined [5f|. To determine the 
confidence region of the parameters of the model, we use 
the following \ 2 function: 
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[i?», ? :- C (r, 4 /r)(W3+3/7)] 



(20) 



where r,i = r(Zi, Q m o = 0.1, Q\o = 0), <7j is the combi- 
nation of the errors in (D) and -D*, and c is a parame- 
ter that minimizes the x 2 for fixed values of the cosmo- 
logical parameters. In our computation, we marginalize 
over P assuming that it is Gaussian distributed such that 
= 1.75 ±0.25 

In Fig. |2| we show contours of constant likelihood 
(95.4% and 68.3%) arising from the radio galaxies test 
for the quartessence Chaplygin model. We can see that 
the FR lib radio galaxies test also does not strongly con- 
strain the QCM models. Only models with A < 0.2 are 
excluded by the data at the 95.4% confidence level. How- 
ever, as we shall see in Sec lVllI when we combine this 
test with the strong lensing test, we get results similar to 
those obtained from SNIa. 



VI. X-RAY GAS MASS FRACTION IN 
GALAXY CLUSTERS 

In recent years, considerable efforts have been devoted 
to determining the matter content of clusters of galax- 
ies. Clusters of galaxies are the most recent large-scale 
structures formed and are also the largest gravitationally 
bound systems known. Therefore, the determination of 
their matter content is quite important because cluster 
properties should approach those of the Universe as a 
whole. A powerful method based on this idea is to mea- 
sure the baryon mass fraction f2{,/f2 m in rich clusters. 
By combining this ratio with Q}, determinations from pri- 
mordial nucleosynthesis, strong constraints on f2 m can be 
placed |5^. As we shall show, this is especially interest- 
ing in quartessence models such as QCM; since the dark 
sector is unified in these models, any strong constraint 
on effective dark matter translates directly into a strong 
constraint on effective dark energy. 

Here we use the method and data of Allen, Schmidt, 
and Fabian 1591 and Allen et al. |60|. These authors 
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In Fig. 4, we show the 68.3% and 95.4% confidence 
contours on the parameters a and A determined from 
the Chandra data. The best-fit value, the solid dot in 
the figure, is located at a ~ and A ~ 0.73 and, as 
discussed in Sec|n] corresponds to the ACDM limit of 
QCM. It is clear from the figure that this test is much 
more restrictive than the others discussed previously. 



FIG. 3: Confidence contours in the (a, A) plane allowed by a 
set of FR lib radio galaxies. 



extracted from Chandra observations the x-ray gas mass 
fraction / gas of nine massive, dynamically relaxed galaxy 
clusters, with rcdshifts in the range 0.08 < z < 0.47, 
and that have converging / gas within a radius ^2500 (ra- 
dius encompassing a region with mean mass density 2500 
times the critical density of the Universe at the cluster 
redshift). ^ 

To determine the confidence region of the parameters 
of the model, we use the following x 2 function in our 
computation: 



[/gas ( z i) 



ft 



gas,! J 



(21) 



7. 



where Zi, / ga s,i; and o"/ gM t arc, respectively, the redshifts, 
the SCDM (h = 0.5) best-fitting values, and the symmet- 
ric root-mean-square errors for the nine clusters as given 
in Ipll and |gjj. In Eq. gTJ, /™° d is the model function 



ig m as d M = 



dl ds 



3/2 



(1 + 0.19>/fe)fi£ \°- 5rf 



a, A 



(22) 



Here, <1a is the angular diameter distance to the cluster, 
fi^f, given by Eq. IjlUfl. is the effective matter density pa- 
rameter, and b is a bias factor that takes into account the 
fact that the baryon fraction in clusters could be lower 
than for the Universe as a whole. In our computations, 
we marginalize over the bias factor assuming that it is 
Gaussian-distributed with b = 0.93 ± 0.05 as suggested 
by gas-dynamical simulations 2 |60t l6ll | . 
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FIG. 4: Constant confidence contours (68.3% and 95.4%) in 
the (a, A) plane determined from the x-ray gas mass fraction 
in nine galaxy clusters from Chandra data. The best-fit value 
is indicated by a dot at the center of the contours. 



[2] We also considered b = 0.824 ± 0.033 for the bias factor (see 



l62ll63l ). In this case, the resulting contours are slightly thinner 
than those in Fig. 4, and the best-fit value for the parameters 
occurs at a ~ -0.12 and A ~ 0.72. 
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FIG. 5: Contours of constant f2„ , corresponding to f2„ — 
0.25 and fi^f = 0.35. 

Independent constraints on Chaplygin models from 
galaxy clusters x-ray data have been presented recently 
in Ref. [64j. In that work, different constraints were ob- 
tained; for instance, the contours are quite insensitive to 
the parameter A. In fact, those authors were interested 
in quintessence Chaplygin, that is, models in which the 
Chaplygin component behaves only as dark energy, while 
here we are considering quartessence Chaplygin (no dark 
matter). Therefore, since Ref. [64| started from different 
priors and did not consider that the Chaplygin compo- 
nent can cluster, it is natural that it reached different 
conclusions. 

It is worth noting that the contours that we have ob- 
tained with the Chandra data correspond to constraints 
on f2°f. To illustrate this, we plot in Fig. 5 the lines 
fi^f = 0.25 and fi^f = 0.35. The similarity of the con- 
tours is evident. We remark that constraints on the QCM 
models similar to those obtained here should be expected 
from large-scale structure data (assuming entropy per- 
turbations as in |l6| ) that are very sensitive to r e ff, the 
QCM effective shape parameter of the mass power spec- 
trum [Hi El, 

r eS = Q^he X p(-Sl h -^^] . (23) 



VII. COMBINED ANALYSIS AND DISCUSSION 

Here we summarize our results, presenting a combined 
analysis of the constraints discussed in the previous sec- 
tions. In Fig. 6, we display the allowed region of the 
parameters A and a from a combination of data from 
FR lib galaxies and gravitational lensing. Although A is 
better constrained than in each one of these experiments 



FIG. 6: Combined results from the statistics of gravitational 
lensing and radio galaxies data. We display the 68.3% and 
95.4% confidence levels in the (a, A) plane. 



separately, it is still not possible to constrain a within our 
expected interval. Even if we include type la supernovae 
data, as in Fig. 7, we are not able to constrain this pa- 
rameter (although the best-fit value is already inside the 
expected interval); only A can be fairly well constrained 
with this set. However, the inclusion of cluster x-ray gas 
mass fraction data, combined with the previous three, as 
in Fig. 8, places significant constraints on both A and 
a. These are the tightest constraints on the QCM pa- 
rameters set to date from the background geometry. The 
allowed interval of a, —0.5 < a < 1, falls within the 
expected interval discussed in Sec. [H] Intriguingly, the 
best-fit value of a is near the ACDM limit of the model, 
a = 0. 

It is seen that QCM models with a > —0.5 arc con- 
sistent with the observables considered here. Moreover, 
even the Chaplygin (a = 1) model — the most theoreti- 
cally motivated of the QCM family — cannot be ruled out 
by current observational data. 

The aim of this paper was to verify whether the 
quartessence Chaplygin model is consistent with cur- 
rent data on the background geometry, and to set con- 
straints on the model parameters. As we have discussed 
in Sec||J even if the ACDM model is in agreement with 
the data, it is still interesting to look for alternative mod- 
els, both from a philosophical point of view and also to 
test the robustness of the model against observational 
data, i.e., to sense to which extent we can modify the 
ACDM paradigm and still be in agreement with the data. 

We argue that, as a model for dark energy (explana- 
tion for accelerated expansion), the Chaplygin gas is not 
particularly attractive. The most promising feature of 
it is in the context of quartessence (unification of dark 
matter and dark energy) . Here we have focused on the re- 
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FIG. 7: Combined results from supernovae, gravitational lens- 
ing statistics, and radio galaxies data. We display the 68.3% 
and 95.4% confidence levels in the (a, ^4) plane. The dot is the 
best-fitting value for the combination of these experiments. 
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FIG. 8: Confidence levels (68.3% and 95.4%) from the com- 
bination of all the observables considered in this paper (su- 
pernovae, gravitational lensing statistics, radio galaxies, and 
cluster x-ray gas mass fraction). The best-fit value is located 
close to the ACDM limit of the model. 



alization of the quartessence scenario with a fluid whose 
background equation of state is given byEq. ©. Rather 
than favoring a specific model, our results show that al- 
ternatives to ACDM are consistent with the data. 

We have addressed only observables that probe essen- 
tially the distance-redshift relation, and it is clear that 
one should look for independent tests, such as, for in- 
stance, in the large-scale structure of the Universe. In 
fact, the QCM has a rich behavior regarding the density 
perturbations. For adiabatic perturbations,the Chaply- 
gin power spectrum has strong oscillations However, 
for positive values of a, the baryon power spectrum is well 
behaved jg^- F° r some entropy perturbations, the QCM 
power spectrum itself is well behaved, even for negative 
values of a [l6|j. Therefore, one should investigate fur- 
ther the large-scale distribution of the dark component 
to better constrain the QCM. In particular, the nonlinear 
regime of structure formation in QCG models should be 
more studied. This should shed some light on the issue of 
the separation into clustered CDM-like and low-density 
negative pressure regions. The impact on the CMB must 
also be more studied [l5t \6ll IfM Ir39l| . 

Of course our analysis could be improved in several as- 
pects. For instance, in the lensing analysis, radio sources 
could be included. Also, we fixed some parameters to 
their best-fitting values, and they could be marginalized 
to provide more robust constraints. However, we expect 
that this will not change qualitatively the results pre- 
sented here. 

The QCM still seems to be a promising model for uni- 
fying dark matter and dark energy. More generically, the 
idea of quartessence has to be explored further, both from 
the particle physics point of view — to search for a first 
principles motivation to it — as well as from the empirical 
side, to constrain quartessence models from observational 
data. 
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